Defence  Research  and  Recherche  et  developpement 
Development  Canada  pour  la  defense  Canada 


Copy  No. 


1*1 


DEFENCE 

Modelling  the  electromagnetic  scattering 
by  3-D  objects  in  a  layered  media 

Electromagnetic  active  detection  of  metallic  targets 


Marius  S.  Birsan 


Defence  RSrD  Canada  -  Atlantic 

Technical  Memorandum 
DRDC  Atlantic  TM  2004-246 
December  2004 


Canada 


Report  Documentation  Page 


Form  Approved 
OMB  No.  0704-0188 


Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 
VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 


1.  REPORT  DATE 

DEC  2004 


2.  REPORT  TYPE 


3.  DATES  COVERED 


4.  TITLE  AND  SUBTITLE 

Modeling  the  electromagnetic  scattering  by  3-D  objects  in  a  layered 
media  (U) 

6.  AUTHOR(S) 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Defence  R&D  Canada  -Atlantic, PO  Box  1012, Dartmouth,  NS,CA,B2Y 
3Z7 


5a.  CONTRACT  NUMBER 


5b.  GRANT  NUMBER 


5c.  PROGRAM  ELEMENT  NUMBER 


5d.  PROJECT  NUMBER 


5e.  TASK  NUMBER 


5f.  WORK  UNIT  NUMBER 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


9.  SPONSORING/MONITORING  AGENCY  NAME(S )  AND  ADDRESS(ES )  10.  SPONSOR/MONITOR' S  ACRONYM(S) 

11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

The  original  document  contains  color  images. 

14.  ABSTRACT 

The  underwater  electromagnetic  (EM)  active  detection  technique  attempts  to  reveal  a  stealth  iron-hull 
vessel  (submarine  or  surface  warship)  located  in  a  horizontal  layer  of  seawater.  Our  goal  is  to  numerically 
investigate  the  EM  scattering  process  in  the  presence  of  a  metallic  vessel  of  arbitrary  shape  placed  inside  a 
lossy  media.  To  solve  the  EM  scattering  problem  involving  perfectly  conducting  (PEC)  objects  we  make  use 
of  the  Method  of  Moments  (MoM)  integral  equation  technique.  The  EM  scattering  numerical  code  was 
written  in  Fortran.  The  solution  is  applied  to  the  investigation  of  EM  scattering  by  a  submarine  modelled 
as  a  tapered  long  cylinder  plus  the  tower  when  is  placed  in  different  positions  relative  to  the  source  and  the 
sensors.  It  is  shown  that,  by  an  appropriate  design  of  the  experiment,  the  signal  indicating  the  presence  of 
the  target  can  be  increased  relative  to  the  signal  obtained  in  its  absence. 


15.  SUBJECT  TERMS 


16.  SECURITY  CLASSIFICATION  OF: 


a.  REPORT 

unclassified 


b.  ABSTRACT 

unclassified 


c.  THIS  PAGE 

unclassified 


17.  LIMITATION  OF 

18.  NUMBER 

ABSTRACT 

OF  PAGES 

38 

19a.  NAME  OF 
RESPONSIBLE  PERSON 


Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


This  page  intentionally  left  blank. 


Modelling  the  electromagnetic  scattering 
by  3-D  objects  in  a  layered  media 

Electromagnetic  active  detection  of  metallic  targets 


Marius  S.  Birsan 


Defence  R&D  Canada  -  Atlantic 

Technical  Memorandum 
DRDC  Atlantic  TM  2004-246 


December  2004 


Author 


Dr.  Ron  Kuwahara 
S/H  SIGNATURES 


Approved  for  release  by 


)RP  Chair 


©  Her  Majesty  the  Queen  as  represented  by  the  Minister  of  National  Defence,  2004 
©  Sa  majeste  la  reine,  representee  par  le  ministre  de  la  Defense  nationale,  2004 


Abstract 


The  underwater  electromagnetic  (EM)  active  detection  technique  attempts  to 
reveal  a  stealth  iron-hull  vessel  (submarine  or  surface  warship)  located  in  a  horizontal 
layer  of  seawater.  It  is  supposed  that  the  vessel  has  no  signature,  acoustic  or 
electromagnetic,  and  the  only  way  to  detect  its  presence  is  by  measuring  the  scattering 
signal  caused  by  an  active  source  placed  also  in  the  water.  Our  goal  is  to  numerically 
investigate  the  EM  scattering  process  in  the  presence  of  a  metallic  vessel  of  arbitrary 
shape  placed  inside  a  lossy  media. 

A  frequently  occurring  geometry  for  the  marine  environment  is  a  horizontal 
stratification  with  constant  conductivity,  permittivity  and  permeability  in  each  layer. 
The  EM  scattering  from  an  arbitrary  shaped  electrical  conducting  object  located  in  a 
layered  medium  is  an  intricate  problem  due  to  the  interaction  between  the  boundaries 
of  the  media  and  the  scattering  object. 

To  solve  the  EM  scattering  problem  involving  perfectly  conducting  (PEC) 
objects  we  make  use  of  the  Method  of  Moments  (MoM)  integral  equation  technique. 
The  scattering  fields  are  expressed  as  surface  integrals  on  the  scattering  object  surface. 
The  integrands  of  these  integrals  consist  of  scalar  products  between  an  unknown 
surface  current  density  and  an  electric  or  magnetic  dyadic  Green's  function  for  the 
electric  and  magnetic  fields,  respectively.  Then  the  unknown  surface  current  density  is 
obtained  from  an  integral  equation  derived  from  the  magnetic  scattered  field 
representation.  The  method  of  moments  (MoM)  is  used  to  transform  the  magnetic 
field  integral  equation  into  a  linear  system  of  equations  solved  for  the  induced  electric 
surface  currents.  Then  the  fields  produced  by  the  equivalent  electric  dipoles 
distributed  on  the  object  surface  are  calculated  by  taking  into  account  the  actual 
layered  structure. 

The  EM  scattering  numerical  code  was  written  in  Fortran.  The  solution  is 
applied  to  the  investigation  of  EM  scattering  by  a  submarine  modelled  as  a  tapered 
long  cylinder  plus  the  tower  when  is  placed  in  different  positions  relative  to  the  source 
and  the  sensors.  It  is  shown  that,  by  an  appropriate  design  of  the  experiment,  the 
signal  indicating  the  presence  of  the  target  can  be  increased  relative  to  the  signal 
obtained  in  its  absence. 
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Resume 


La  technique  active  de  detection  electromagnetique  sous-marine  vise  la 
decouverte  des  navires  furtifs  a  coque  d’acier  (sous-marin  ou  vaisseau  de  surface), 
places  dans  une  couche  horizontale  d’eau  de  mer.  On  presume  que  le  navire  n’a  pas  de 
signature  acoustique  ou  electromagnetique  et  que  la  seule  fagon  de  deceler  sa  presence 
est  de  mesurer  la  diffusion  d’un  signal  emis  par  une  source  immergee.  Notre  objectif 
etait  d’etudier  numeriquement  la  diffusion  electromagnetique  produite  par  un  navire 
metallique  de  forme  arbitraire  place  dans  un  milieu  a  perte  et  stratifie. 

Generalement,  on  peut  representer  l’environnement  marin  comme  une  serie  de 
strates  horizontales,  caracterisees  par  une  conductivity,  une  permittivity  et  une 
permeabilite  constantes.  La  diffusion  electromagnetique  d’un  objet  conducteur  de 
forme  arbitraire  place  dans  un  milieu  stratifie  est  un  probleme  difficile,  etant  donne 
l’interaction  entre  les  limites  du  milieu  et  l’objet  diffusant. 

Nous  avons  resolu  le  probleme  de  la  diffusion  electromagnetique  par  des 
objets  parfaitement  conducteurs  en  exploitant  la  methode  des  moments  pour  resoudre 
1’ equation  integrale.  Les  champs  diffuses  sont  exprimes  par  des  integrates  de  surface, 
a  la  surface  des  objets  diffusants.  Les  fonctions  a  integrer  sont  les  produits  scalaires 
d’un  courant  de  surface  inconnu  par  une  fonction  dyadique  de  Green,  electrique  ou 
magnetique,  selon  le  champ  considere.  Ensuite,  nous  obtenons  la  densite  de  courant 
inconnue,  a  partir  de  l’equation  integrale  derivee  de  la  representation  du  champ 
magnetique  diffuse.  L’appbcation  de  la  methode  des  moments  pennet  de  transformer 
1’ equation  integrale  du  champ  magnetique  en  un  systeme  d’ equations  lineaires  que 
Ton  resout  pour  trouver  les  courants  de  surface  induits.  Nous  calculons  alors  les 
champs  produits  par  les  dipoles  electriques  equivalents  distribues  a  la  surface  de 
l’objet,  en  tenant  compte  de  la  structure  stratifiee. 

Nous  avons  ecrit  le  programme  de  simulation  de  la  diffusion 
electromagnetique  en  FORTRAN.  Nous  avons  utilise  la  solution  pour  etudier  la 
diffusion  electromagnetique  par  un  sous-marin  modelise  comme  un  long  cylindre 
effile  surmonte  d’un  kiosque,  place  dans  differentes  positions  relativement  a  la  source 
et  aux  capteurs.  Nous  montrons  grace  a  une  conception  appropriee  de  l’experience, 
que  le  signal  indiquant  la  presence  d’une  cible  peut  etre  superieur  au  signal  regu  en 
son  absence. 
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Executive  summary 


Introduction 

For  several  years,  the  active  electromagnetic  (EM)  detection  method  has  been 
experimented  in  our  laboratory  as  a  possible  surveillance  system  for  surface  and 
submerged  vessels.  The  experiments  were  motivated  by  the  potential  of  the 
electromagnetic  active  detection  systems  to  outperform  the  acoustic  systems  in  certain 
situations.  The  modem  surface  and  submerged  vessels  are  very  quiet  acoustically,  so 
that  they  are  difficult  to  detect  by  passive  acoustic  methods.  Moreover,  in  shallow 
waters  and  surf  zones  both  passive  and  active  acoustic  detection  techniques  perform 
poorly  due  to  multiple  scattering  at  the  boundaries  and  high  noise  level. 

The  underwater  EM  active  detection  technique  attempts  to  reveal  a  stealth 
iron-hull  vessel  (submarine  or  surface  warship)  located  in  a  horizontal  layer  of 
seawater.  It  is  supposed  that  the  vessel  has  no  acoustic  or  EM  signature  and  the  only 
way  to  detect  its  presence  is  by  measuring  the  scattering  signal  caused  by  an  active 
source  placed  also  in  the  water.  The  conducting  water  layer  is  situated  between  a  half¬ 
space  of  air  and  a  known  number  of  layers  representing  sediments  of  lower  electrical 
conductivity.  The  electrical  conductivity  of  the  metallic  hull  is  at  least  seven  orders  of 
magnitude  higher  than  the  conductivity  of  the  surrounding  media,  so  that  the  vessel 
could  be  considered  as  a  perfectly  conducting  3-D  object.  An  electric  or  magnetic 
source  placed  in  the  water  creates  an  EM  field  that  propagates  in  the  layered  marine 
environment  and  illuminates  the  vessel,  which  in  turn  will  generate  a  scattered  signal. 

The  ability  to  set  up  appropriate  EM  experiments  and  to  enhance  the 
interpretation  of  the  measured  signals  makes  it  necessary  to  use  realistic  and  efficient 
numerical  models.  Our  goal  is  to  numerically  investigate  the  scattering  process  in  the 
presence  of  the  metallic  vessel.  This  report  presents  a  mathematical  model  of  the  EM 
scattering  caused  by  a  perfectly  conducting  (PEC)  body  of  arbitrary  shape  placed 
inside  a  lossy-layered  media. 

Results 

To  solve  the  EM  scattering  problem  involving  perfectly  conducting  (PEC) 
objects  in  layered  media  we  make  use  of  the  Method  of  Moments  (MoM)  integral 
equation  technique.  The  scattering  fields  are  expressed  as  surface  integrals  on  the 
scattering  object  surface.  The  integrands  of  these  integrals  consist  of  scalar  products 
between  an  unknown  surface  current  density  and  an  electric  or  magnetic  dyadic 
Green's  function  for  the  electric  and  magnetic  fields,  respectively.  The  dyadic  Green's 
functions  are  fundamental  wave  propagation  solutions  to  a  certain  stratified  geometry. 
In  other  words,  the  interaction  between  the  scattering  object  and  the  air-seawater  and 
seawater-seabed  interfaces  is  taken  into  account.  Then  the  unknown  surface  current 
density  is  obtained  from  an  integral  equation  of  the  second  kind  derived  from  the 
magnetic  scattered  field  representation.  This  approach  for  solving  the  scattering 
problem  is  called  the  magnetic  field  integral  equation  (MFIE).  The  method  of 
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moments  (MoM)  is  used  to  transform  the  magnetic  field  integral  equation  (MFIE)  into 
a  linear  system  of  equations.  The  resulting  matrix  equation  is  solved  for  the  induced 
electric  surface  currents  by  Gauss  elimination.  Then  the  fields  produced  by  the 
equivalent  electric  dipoles  distributed  on  the  object  surface  are  calculated  by  taking 
into  account  the  actual  layered  structure.  The  EM  scattering  numerical  code  was 
written  in  Fortran.  The  solution  is  applied  to  investigate  the  EM  scattering  by  a 
submerged  submarine  when  it  is  placed  in  different  positions  relative  to  the  source 
and  the  sensors.  The  measured  electric  and  magnetic  fields  are  composed  as  the  sum 
of  the  direct  (primary)  incoming  fields  from  the  source,  i.e.  the  wave  propagation 
fields  without  any  scattering  object  in  the  environment,  and  the  scattered  (secondary) 
fields  due  to  the  object.  It  is  shown  that,  by  an  appropriate  design  of  the  experiment, 
the  signal  indicating  the  presence  of  the  target  can  be  increased  relative  to  the  signal 
obtained  in  its  absence. 

Significance 

The  mathematical  model  of  the  EM  scattering  by  metallic  objects  submerged  in 
seawater  or  buried  in  the  sediment  offers  useful  information  for  designing  the 
underwater  EM  active  detection  (of  vessels  or  mines)  experiments  and  interpreting  the 
experimental  data. 

Future  plans 

In  the  active  detection  of  the  vessels,  the  EM  source  and  sensor  are  located  in 
the  seawater,  but  a  particular  advantage  of  inductive  coupling  is  that  it  permits  the  use 
of  EM  systems  in  aircraft  as  an  additional  tool  to  the  Magnetic  Anomaly  Detection 
(MAD)  system. 

Because  the  direct  (primary)  field  is  much  higher  than  the  secondary  field,  it  is 
necessary  to  use  measurement  techniques  to  diminish  the  primary  field  effect.  Such 
methods  were  developed  and  electromagnetic  equipment  has  been  widely  applied  in 
mineral  exploration  reconnaissance.  For  the  future,  it  is  intended  to  verify  the 
applicability  of  the  geophysical  electromagnetic  methods  in  underwater  active 
detection  and  to  develop  new  measurement  techniques,  in  the  water  or  in  the  air,  if 
necessary. 


Birsan  M.  2004.  Modelling  the  electromagnetic  scattering  by  3-D  objects  in  a  layered 
media.  DRDC  Atlantic  TM  2004-246.  DRDC  Atlantic. 
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Sommaire 


Introduction 

Depuis  plusieurs  annees,  nous  etudions  en  laboratoire  la  possibilite  d’utiliser  la 
detection  electromagnetique  active  pour  surveiller  les  navires  immerges  ou  en  surface. 
Ces  experiences  sont  motivees  par  la  possibilite  que,  dans  certaines  situations,  les 
systemes  de  detection  electromagnetique  active  soient  superieurs  aux  systemes 
acoustiques.  Etant  donne  la  grande  discretion  acoustique  des  navires  modemes  —  tant 
en  surface  ou  sous  l’eau  —  il  est  devenu  difficile  de  les  detecter  par  les  methodes 
acoustiques  passives.  En  outre,  dans  les  hauts  fonds  et  les  zones  de  deferlement,  les 
techniques  acoustiques  de  detection,  actives  ou  passives,  sont  inefficaces,  a  cause  des 
diffusions  multiples  sur  les  bords  et  l’intensite  du  bruit  ambiant. 

La  technique  active  de  detection  electromagnetique  sous-marine  vise  la 
decouverte  des  navires  furtifs  a  coque  d’acier  (sous-marin  ou  vaisseau  de  surface), 
places  dans  une  couche  horizontale  d’eau  de  mer.  On  presume  que  le  navire  n’a  pas  de 
signature  acoustique  ou  electromagnetique  et  que  la  seule  fag  on  de  deceler  sa  presence 
est  de  mesurer  la  diffusion  d’un  signal  emis  par  une  source  immergee.  La  couche 
d’eau  conductrice  s’etend  entre  un  demi-espace  d’air  et  le  lit  forme  d’un  nombre 
connu  de  strates  de  sediments,  dont  la  conductivite  electrique  est  plus  basse.  La 
conductivite  electrique  d’une  coque  metallique  est,  au  moins,  sept  ordres  de 
magnitude  plus  elevee  que  celle  du  milieu  ambiant  et,  ainsi,  on  peut  considerer  que  le 
navire  est  un  objet  tridimensionnel  parfaitement  conducteur.  Une  source  electrique  ou 
magnetique  placee  dans  l’eau  genere  un  champ  electromagnetique  qui  se  propage  dans 
le  milieu  marin  stratifie  et  illumine  le  vaisseau  qui  diffuse  un  signal. 

Puisque  nous  pouvons  realiser  des  experiences  electromagnetiques  adequates  et 
accentuer  les  signaux  mesures  pour  les  interpreter,  nous  devons  utiliser  des  modeles 
numeriques  realistes  et  efficaces.  Notre  objectif  etait  done  d’etudier  numeriquement  la 
diffusion  par  un  navire  metallique.  Dans  ce  rapport,  nous  presentons  un  rnodele 
mathematique  de  la  diffusion  electromagnetique  par  un  corps  parfaitement 
conducteur,  de  forme  arbitraire,  place  dans  un  milieu  a  perte  et  stratifie. 

Resultats 

Nous  avons  resolu  le  probleme  de  la  diffusion  electromagnetique  par  des 
objets  parfaitement  conducteurs  dans  un  milieu  stratifie,  en  exploitant  la  methode  des 
moments  pour  resoudre  l’equation  integrale.  Les  champs  diffuses  sont  exprimes  par 
des  integrales  de  surface,  a  la  surface  des  objets  diffusants.  Les  fonctions  a  integrer 
sont  les  produits  scalaires  d’un  courant  de  surface  inconnu  par  une  fonction  dyadique 
de  Green,  electrique  ou  magnetique,  selon  le  champ  considere.  Les  fonctions 
dyadiques  de  Green  sont  les  solutions  fondamentales  de  propagation  dans  une 
geometrie  stratifiee  donnee.  En  d’autres  tennes,  il  faut  considerer  l’interaction  entre 
1’ objet  diffusant  et  les  interfaces  air-eau  et  eau-fond.  Ensuite,  nous  ob tenons  la  densite 
de  courant  inconnue,  a  partir  de  l’equation  integrale  de  deuxieme  espece,  derivee  de  la 


VI 


DRDC  Atlantic  TM  2004-246 


representation  du  champ  magnetique  diffuse.  On  appelle  equation  integrale  du  champ 
magnetique,  cette  methode  de  resolution  du  probleme  de  diffusion.  L’application  de  la 
methode  des  moments  nous  permet  de  transformer  1’ equation  integrale  du  champ 
magnetique  en  un  systeme  d’ equations  lineaires.  Nous  utilisons  ensuite  la  methode 
d’ elimination  de  Gauss  pour  resoudre  1’ equation  matricielle  resultante  pour  les 
courants  de  surface  induits. 

Nous  calculons  alors  les  champs  produits  par  les  dipoles  electriques 
equivalents  distribues  a  la  surface  de  l’objet,  en  tenant  compte  de  la  structure 
strati  lice.  Nous  avons  ecrit  le  programme  de  simulation  de  la  diffusion 
electromagnetique  en  FORTRAN.  La  solution  est  utilisee  pour  etudier  la  diffusion 
electromagnetique  par  un  sous-marin  immerge,  place  dans  differentes  positions 
relativement  a  la  source  et  aux  capteurs.  Les  champs  electriques  et  magnetiques 
mesures  sont  la  resultante  des  champs  directs  (primaires)  emis  par  la  source  —  soit  les 
champs  de  propagation  des  ondes  dans  un  environnement  ne  comportant  pas  d’objet 
diffusant  —  et  des  champs  diffuses  par  l’objet  (secondaires).  Nous  montrons  grace  a 
une  conception  appropriee  de  l’experience,  que  le  signal  indiquant  la  presence  d’une 
cible  peut  etre  superieur  au  signal  rcgu  en  son  absence. 

Portee 

Le  modele  mathematique  de  diffusion  electromagnetique  par  des  objets 
metalliques  plonges  dans  l’eau  de  mer  ou  enfouis  dans  les  sediments  fournit  des 
informations  utiles  pour  la  conception  des  experiences  de  detection  electromagnetique 
active  sous-marine  (de  navires  ou  de  mines)  et  pour  interpreter  les  donnees 
experimentales. 

Futures  recherches 

Pour  la  detection  active  des  navires,  la  source  electromagnetique  et  les  capteurs 
sont  situes  dans  l’eau  de  mer.  Toutefois,  le  couplage  inductif  presente  l’avantage 
supplemental  de  permettre  1 ’utilisation  de  systemes  electromagnetiques  aeroportes 
comme  outil  additionnel  au  systeme  de  detection  d’anomalie  magnetique. 

Puisque  le  champ  direct  (primaire)  est  beaucoup  plus  intense  que  le  champ 
secondaire,  il  est  necessaire  d’utiliser  des  techniques  de  mesure  pennettant  de  reduire 
l’effet  du  champ  primaire.  Ces  methodes  ont  deja  ete  mises  au  point  et  on  a  largement 
utilise  l’equipement  electromagnetique  issu  de  ces  travaux  pour  la  prospection 
mini  ere.  A  l’avenir,  nous  prevoyons  experimenter  des  methodes  electromagnetiques 
geophysiques  pour  la  detection  sous-marine  active  et  la  mise  au  point  de  nouvelles 
techniques  de  mesure,  sous  l’eau  ou  dans  Pair,  au  besoin. 


Birsan  M.  2004.  Modelling  the  electromagnetic  scattering  by  3-D  objects  in  a  layered 
media.  (Modelisation  de  la  diffusion  electromagnetique  par  des  objets  tridimensionnels 
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1.  Introduction 


The  detection  and  classification  of  submerged  objects  in  seawater  is  mainly 
performed  using  hydro-acoustic  or  electromagnetic  methods.  These  methods  can 
further  be  divided  into  passive  and  active  methods.  A  passive  system  consists  of 
sensors  that  detect  either  the  sounds  emitted  by  the  submerged  objects  or  their 
electromagnetic  signature.  An  active  system  consists  of  both  sources  and  sensors.  The 
sources  radiate  acoustic  or  electromagnetic  energy  towards  the  object  and  the  sensors 
are  designed  to  detect  any  reflected  energy. 

For  several  years,  the  active  electromagnetic  (EM)  detection  method  has  been 
experimented  in  our  laboratory  as  a  possible  surveillance  system.  The  experiments 
were  motivated  by  the  potential  of  the  electromagnetic  active  detection  systems  to 
outperform  the  acoustic  systems  in  certain  situations.  The  modem  surface  and 
submerged  vessels  are  very  quiet  acoustically,  so  that  they  are  difficult  to  be  detected 
by  passive  acoustic  methods.  Moreover,  in  shallow  waters  and  surf  zones  both  passive 
and  active  acoustic  detection  techniques  perform  poorly  due  to  multiple  scattering  at 
the  boundaries  and  high  noise  level.  It  is  also  assumed  that,  in  the  future,  the  vessels 
will  be  equipped  with  degaussing  and  anti-corrosion  systems  designed  to  completely 
cancel  their  magnetic  and  electric  signatures. 

The  active  electromagnetic  (EM)  detection  of  the  surface  (partially 
submerged)  and  submerged  vessels  in  the  seawater  or  of  the  buried/partially  buried 
objects  in  the  sediment  is  based  on  the  measurement  of  the  response  (scattered)  signal 
from  the  metallic  object  due  to  an  EM  source.  The  method  involves  the  propagation  of 
time-varying,  low  frequency  electromagnetic  fields  in  and  over  the  ocean.  The 
detector  receives  its  signal  by  induction.  The  measured  total  electric  and  magnetic 
fields  are  composed  as  the  sum  of  the  direct  incoming  fields  from  the  source,  i.e.  the 
wave  propagation  fields  without  any  scattering  object  in  the  environment,  and  the 
scattered  fields  due  to  the  object.  In  the  active  detection  of  the  vessels,  the  EM  source 
and  sensor  are  located  in  the  seawater,  but  a  particular  advantage  of  the  inductive 
coupling  is  that  it  permits  the  use  of  EM  systems  in  aircraft.  Airborne  electromagnetic 
in  combination  with  magnetic  equipment  has  been  widely  applied  in  mineral 
exploration  reconnaissance,  so  that  it  could  be  an  additional  tool  to  the  Magnetic 
Anomaly  Detection  (MAD)  system. 

The  ability  to  set  up  appropriate  EM  experiments  and  to  enhance  the 
interpretation  of  the  measured  signals  makes  it  necessary  to  use  realistic  and  efficient 
numerical  models.  It  is  the  purpose  of  this  report  to  present  a  mathematical  model  of 
the  EM  scattering  caused  by  a  perfectly  conducting  (PEC)  body  of  arbitrary  shape 
placed  inside  a  lossy-layered  media. 

To  solve  the  EM  scattering  problem  involving  perfectly  conducting  (PEC) 
objects  we  make  use  of  the  Method  of  Moments  (MoM)  integral  equation  technique. 
The  methodology  used  in  this  report  follows  closely  the  one  used  in  the  well-known 
code  NEC  (Numerical  Electromagnetic  Code)  developed  at  the  Lawrence  Livennore 
Laboratory  for  the  solution  of  the  EM  scattering  and  radiation  problems  involving 
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objects  of  arbitrary  shape  in  free  space.  The  NEC  Fortran  code  is  available  in  the 
public  domain.  However,  for  the  underwater  active  detection  applications,  it  is 
necessary  to  formulate  the  scattering  problem  for  an  arbitrary  shaped  PEC  object 
located  in  a  layered  media,  which  is  an  intricate  scattering  problem.  The  difficulty  of 
such  a  problem  is  due  to  the  interaction  between  the  boundaries  of  the  media  and  the 
scattering  object. 

A  frequently  occurring  geometry  for  the  marine  environment  is  a  horizontal 
stratification  with  constant  conductivity,  permittivity  and  permeability  in  each  layer. 
Examples  of  papers  dealing  with  electromagnetic  very  low  frequency  (VLF) 
scattering  and  wave  propagation  in  such  environmental  configurations  are  references 
[1-6].  The  mathematical  techniques  used  for  the  scattering  part  of  the  problem  are 
mostly  various  integral  equation  (IE)  fonnulations  with  different  method  of  moments 
(MoM)  concepts,  [2,  3,  4,  6].  There  are  also  some  hybrid  techniques  where  the 
integral  equation  and  the  finite  elements  (FE)  method  are  combined,  [7,  8],  The  FE 
method  is  used  to  model  an  inhomogeneous  scattering  object  and  the  integral  equation 
technique  for  taking  care  of  the  surroundings  to  the  object.  The  two  methods  are 
linked  together  through  the  boundary  conditions  at  a  circumscribing  surface  to  the 
scattering  object.  Other  methods  that  are  useful  for  modelling  inhomogeneous  objects 
and  non-plane  stratified  environments  are  the  finite  difference  time  domain  (FDTD) 
method  [1],  and  FE  methods  for  both  the  scattering  object  and  the  surrounding 
environment  [9].  An  efficient  method  for  modelling  the  scattering  by  a  heterogeneous 
collection  of  spatially  separated  objects  is  a  modified  recursive  T-matrix  algorithm 
given  in  [10]. 

The  underwater  EM  active  detection  technique  attempts  to  reveal  a  stealth 
iron-hull  vessel  (submarine  or  surface  warship)  located  in  a  horizontal  layer  of 
seawater.  It  is  supposed  that  the  vessel  has  no  EM  signature  and  the  only  way  to 
detect  its  presence  is  by  measuring  the  scattering  signal  caused  by  an  active  source 
placed  also  in  the  water.  The  conducting  water  layer  is  situated  between  a  half-space 
of  air  and  a  known  number  of  layers  representing  sediments  of  lower  electrical 
conductivity.  The  electrical  conductivity  of  the  metallic  hull  is  about  seven  orders  of 
magnitude  higher  than  the  conductivity  of  the  surrounding  media,  so  that  the  vessel 
could  be  considered  as  a  perfect  conducting  3-D  object.  The  paint  or  other  insulating 
materials  on  the  hull  do  not  affect  the  perfectly  conducting  (PEC)  assumption  because 
the  wavelength  of  the  EM  signal  is  much  bigger  than  the  thickness  of  the  insulation. 
An  electric  or  magnetic  source  placed  in  the  water  creates  an  EM  field  that  propagates 
in  the  layered  marine  environment  and  illuminates  the  vessel,  which  in  turn  will 
generate  a  scattered  signal.  Our  goal  is  to  numerically  investigate  the  scattering 
process  from  the  vessel.  The  infonnation  is  useful  in  designing  the  EM  active 
detection  experiment  and  interpreting  the  experimental  data. 

Pure  EM-wave  propagation  in  horizontally  stratified  conducting  media  without 
object  scattering  is  treated  in  [11],  The  sources  are  elementary  vertical  and  horizontal 
current  or  magnetic  dipoles  placed  in  an  arbitrary  layer  of  a  N-layers  configuration. 
The  quasi-static  approximation  given  in  [12]  is  used. 
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The  scattering  object  is  a  smooth  perfectly  conducting  3-D  object,  which  can 
be  placed  in  any  of  the  layers  of  an  arbitrary  horizontally  stratified  dissipative 
environment.  The  surface  of  the  object  is  modelled  by  perfectly  conducting  flat 
surface  patches. 

The  scattering  fields  are  expressed  as  surface  integrals  on  the  scattering  object 
surface.  The  integrands  of  these  integrals  consist  of  scalar  products  between  an 
unknown  surface  current  density  and  an  electric  or  magnetic  dyadic  Green's  function 
for  the  electric  and  magnetic  fields,  respectively.  The  dyadic  Green's  functions  are 
fundamental  wave  propagation  solutions  to  a  certain  stratified  geometry.  In  other 
words,  the  interaction  between  the  scattering  object  and  the  air-seawater  and  seawater- 
seabed  interfaces  is  taken  into  account.  The  N-layers  code  is  utilised  to  construct  the 
dyadic  Green’s  functions  when  the  geometry  of  the  environment  is  horizontally 
stratified  in  more  than  two  conducting  layers.  Then  the  unknown  surface  current 
density  is  obtained  from  an  integral  equation  of  the  second  kind  derived  from  the 
magnetic  scattered  field  representation.  This  approach  for  solving  the  scattering 
problem  is  called  the  magnetic  field  integral  equation  (MFIE). 

The  integral  equation  approach  is  best  suited  to  structures  with  dimensions  up 
to  several  wavelengths.  When  the  structure  size  is  increased  relative  to  wavelength, 
standard  high-frequency  approximations  such  as  physical  optics  or  geometric  theory 
of  diffraction  may  be  more  suitable  approaches.  For  the  underwater  active  detection 
applications,  the  source  frequency  is  low  enough,  so  that  the  wavelength  is  usually 
bigger  than  the  largest  dimension  of  the  object  to  be  detected. 

The  method  of  moments  (MoM)  is  used  to  transform  the  magnetic  field 
integral  equation  (MFIE)  into  a  linear  system  of  equations.  The  resulting  matrix 
equation  is  solved  for  the  induced  electric  surface  currents  by  Gauss  elimination.  Next 
the  fields  produced  by  the  equivalent  electric  dipoles  distributed  on  the  object  surface 
are  calculated  by  taking  into  account  the  actual  layered  structure.  The  scattering 
numerical  code  was  written  in  Fortran.  The  solution  is  applied  to  the  investigation  of 
EM  scattering  by  a  submarine  modelled  as  a  tapered  long  cylinder  plus  the  tower 
when  is  placed  in  different  positions  relative  to  the  source  and  the  sensors. 
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2.  The  formulation  of  the  problem 


Consider  a  smooth  perfect  conducting  (PEC)  object  submerged  in  a  horizontal 
layer  of  seawater  of  depth  dw  (Figure  1).  The  X  and  Y  directions  denote  the 
horizontal  plane.  This  is  the  plane  in  which  the  structure  is  unifonn  in  its 
electromagnetic  parameters.  The  Z  direction  is  the  vertical  direction  pointing 
downwards  in  which  the  structure  varies  in  its  properties.  The  origin  of  this  Cartesian 
reference  frame  is  located  on  the  interface  of  air  and  seawater.  A  layered  conducting 
medium  occupies  the  half-space  Z  >  0  of  the  rectangular  coordinate  system  (X,  Y,  Z) 
and  has  the  conductivity  cii  for  0  <  Z  <  dw,  02  for  dw  <  Z  <  dw+ds,  where  ds  is  the 
sediment  thickness,  and  so  on.  The  region  Z  <  0  is  taken  to  be  free-space.  The  position 
of  the  current  source  has  the  coordinates  X  =  Y  =  0,  Z  =  hd  (0  <  hd  <  dw).  The  forcing 
current  has  a  sinusoidal  temporal  dependence,  so  that  the  time  dependence  factor 
exp(icot)  is  suppressed  throughout. 


Figure  1.  Geometry  of  the  EM  scattering  problem  in  a  marine  environment. 
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The  wave  number  k  is  given  by  k  =  ^cojult(co£  +  i a)  ,  where  co  is  the  angular 

frequency,  s  the  pennittivity  and  a  the  conductivity,  and  it  is  assumed  that  all  layers 
are  non-magnetic,  i.e.  the  permeability  is  //0  everywhere.  The  much  higher 
conductivity  of  seawater  compare  to  air  or  sediment  leads  to  a  much  greater  complex 
wave  number  for  EM  waves  in  seawater.  At  low  frequencies,  the  wave  number  in 
water  is  approximated  as: 

kx  =  /?[  +iax  ~  ( icoJu0cTl)u 2  =  (\  +  i\co/u0(7x  12 )1/2  (1) 

The  associated  wavelength  k\  =  2tx/P  i  is  much  shorter  than  that  in  the  air  at  the 
same  frequency.  It  follows  that  a  much  lower  frequency  must  be  used  in  order  to 
formulate  the  scattering  problem  in  tenns  of  integral  equations  valid  when  the 
characteristic  geometric  dimension  of  the  object  is  smaller  than  the  operating 
wavelength.  For  example,  if  the  scattering  object  is  a  submarine  with  a  length  of  50m, 
the  frequency  must  be  smaller  than  225  Hz  (k\  =  105m  for  seawater  with  ai  =  4S/m). 

The  wave  equations  for  the  electric  and  magnetic  fields  Ej  and  Hj,  in  layer  j  (j 
=  0,  1,2,  3),  are  derived  from  the  Maxwell’s  equations  for  time-harmonic  fields: 

V  x  E  =  icou, .  H 

(2) 

VxH,  =Oj  Ej  +J.+VxM. 

where  Jj  and  V  x  Mj  represent  the  external  electric  or  magnetic  source,  and  the 
displacement  currents  were  neglected.  The  boundary  conditions  require  that  the 
magnetic  field  and  the  tangential  component  of  the  electric  field  are  continuous  at 
each  interface  (j  =  0,  1,  2): 

zx(E, -E^,)=0,  ix(H,-H.J=0  (3) 

Pure  EM-wave  propagation  in  horizontally  stratified  conducting  media  without 
object  scattering  is  treated  in  [11].  The  sources  are  elementary  vertical  and  horizontal 
current  or  magnetic  dipoles  placed  within  one  layer  of  a  N-layers  configuration.  The 
quasi-static  approximation  given  in  [12]  is  used. 

The  total  electric  E-  and  magnetic  //-fields  at  the  observation  point  are  the 
sums  of  the  incident  (primary)  fields  from  the  source  and  the  scattered  (secondary) 
fields  from  the  PEC  3-D  object,  that  is  E  =  E'+E'and  H  =  H'+HS.  The 
incident E‘-  and  //'-fields  are  the  wave  propagation  solutions  satisfying  the 
prescribed  horizontally  stratified  geometry  without  the  scattering  object.  The 
scattering  contributions  from  the  object,  /"and  H  ,  also  satisfy  the  N-layers 
geometry.  This  means  that  all  multiple  scattering  effects  are  taken  into  account. 
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In  the  integral  equation  formulation,  the  scattered  fields  ZT  and  Hs  are 
expressed  as  integral  representations  over  the  scattering  object  surface.  For  a  perfect 
conducting  object  the  representations  are: 

E 5  (r)  =  -i  o)ju0  J  G  e  (r,  r)  •  Js  (r')dS(r') 

f  S  (4) 

H'(r)  =  -jGm(r,r')-Js(r')^(r') 

s 

where  Js(r')  is  the  surface  current  density  producing  the  scattering  fields.  The 
electric  and  magnetic  dyadic  Green’s  functions  Ge(r,r')  and  G m (r.r')arc  3x3 
matrices  where  the  columns  are  the  E(r)-  and  //(r)  -fields,  respectively,  from 
nonnalised  electric  current  dipoles  in  three  mutually  perpendicular  directions  at  the 
point  r'in  the  selected  stratified  geometry.  This  means  that  the  expressions  in  (4)  are 
also  solutions  to  that  geometry.  The  vectors  r  and  r'  may  position  the  observation 
point  and  the  source  in  the  same  layer  or  in  different  layers. 

If  the  surface  current  Js(r')  is  induced  by  an  external  incident  field,  H  ,  then 
the  total  magnetic  field  just  inside  the  PEC  object  surface  is  zero: 

H'(r)  +  H'(r)  =  0  (5) 

for  r  just  inside  the  surface  of  the  PEC  object. 

The  unknown  surface  current  density  J^r)  on  the  scattering  object  surface 

may  be  calculated  from  an  integral  equation  of  the  second  kind.  This  integral  equation 
is  obtained  by  operating  with  the  outward  unit  nonnal  n(r)  from  the  left  in  equation 

(5)  and  substituting  the  Hs  -field  representation  from  equation  (4).  Letting  r  -A  S 
from  inside  the  scattering  object  along  the  normal,  the  surface  component  of  equation 
(4)  is: 


- |Js(r)  +  {n(r) x(Gm(r,r  )-Js(r'))^(r')  =  -n(r)xH!(r)  ,  r  eS  (6) 

4  5 

provided  that  the  surface  S  is  smooth. 

The  dyadic  Gm  contains  the  three  magnetic  field  solutions  to  three  volume 
sources  that  are  nonnalized  infinitesimal  current  dipoles  in  three  mutual  perpendicular 
directions  located  at  r’.  The  wave  equation  for  the  dyadic  Green’s  function  with  the 
source  placed  in  water  (j  =  1)  is: 

v X V x G „ (r, r’)- k)  G„. (r,r’|=Vx [i S(r - r')] <5I; 
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where  I  is  the  unit  diagonal  dyadic  and  8(r-r’)  is  the  delta  Dirac  function.  The 
boundary  conditions  at  the  region  interfaces  for  Gm  are  the  same  as  for  the  H  field. 

The  vector  Js(r')  is  tangential  to  the  surface  S.  It  can  therefore  be 

decomposed  as  Js(r')  =  e,(r')J1(r')  +  e2(r')J2(r')  where  ej(r')  and  e2(r')  are  two 
orthogonal  unit  tangential  vectors  at  the  point  r  on  S,  and  (r’)  x  e2 (r’)  =  n(r’) .  The 
vector  integral  equation  (6)  is  then  projected  onto  these  tangential  vectors  at  the  point 
r’  to  obtain  two  scalar  equations.  By  using  the  identity  u  •  (v  x  w)  =  (u  x  v)-  w  and 
noting  that  ej(r’)x  n(r’)  =  -e2(r’)  and  e2(rjxn(rj  =  e1(r'),  the  scalar  equations 
are: 


-  yei  (r) '  J5  (r)  -  J  e2  (r)  •  (Gb  (r,  r)  •  Js  (r'j)dS(r)  =  e2  (r)  •  H  (r) 

(7) 

-e2  (r)  •  J5  (r)  -  J  ei  (r)  •  (G„,  (r,  r)  •  Js  (r'))dS(r)  =  ej  (r)  •  H''  (r) 

These  two  equations  are  sufficient  because  there  is  no  nonnal  component  of  the 
surface  current. 

The  system  of  coupled  scalar  equations  (7)  is  solved  by  the  Method  of 
Moments  approach.  This  means  that  the  surface  is  discretized  and  that  the  unknown 
components  Ji  and  J2  are  expanded  in  a  suitable  set  of  basis  functions.  In  this  way,  the 
integral  equations  are  converted  into  a  system  of  linear  equations  where  the  expansion 
coefficients  are  the  unknown. 
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3.  The  method  of  moments  (MoM) 


Equations  (7)  have  the  form:  L(g)  =  h  ,  where  L  is  a  known  linear  operator,  h 
is  a  known  excitation  function,  and  g  is  the  response  function.  For  equations  (7),  L  is 
an  integral  operator.  The  objective  here  is  to  determine  g  when  L  and  h  are  specified. 
While  the  problem  is  often  intractable  in  closed  form,  the  linearity  of  the  operator  L 
makes  a  numerical  solution  possible.  Such  a  solution  is  provided  by  the  Method  of 
Moments  (MoM)  where  the  unknown  response,  g,  is  represented  as  a  linear 
combination  of  known  basis  or  expansion  functions,  gm(u): 

M 

g{u)=Yja'nSm(U)  (8) 

m= 1 


The  linearity  of  the  operator  L  makes  it  possible  to  write: 

M 

X  L(g,n  )  =  h 

m= 1 


(9) 


For  finite  M  this  equality  is  usually  approximate.  The  choice  of  the  basis  functions  is 
important  for  the  efficiency  and  the  accuracy  of  the  solution.  They  are  selected  so  that 
each  L(gn)  can  be  evaluated  conveniently,  preferably  in  closed  form. 

Equation  (9)  has  M  unknowns.  To  solve  for  the  M  unknowns,  a  set  of  M 
equations  are  obtained  by  introducing  the  testing  or  weighting  functions,  wn,  and 
taking  the  inner  product  (an  integration  over  the  structure  surface  in  this  case)  of  (9): 

M 

Y,0C’n{Wn^L{gm))  =  {Wn^h\  171  =  1,2,3, ...,M  (10) 

m= 1 


The  weighting  functions  must  be  linearly  independent,  so  that  the  M  equations  (10) 
will  be  linearly  independent.  Also  their  selection  is  based  on  the  minimization  of  the 
computational  requirements  to  evaluate  the  inner  product. 

The  system  of  equations  (10)  can  be  written  in  matrix  form: 


[G\[A]  =  [H] 

Gm„  =  ( W,n  >  L(g„  ));  An  =  am  ;  Hn  =  (W n  > ll) 


(ID 


that  could  be  solved  for  [A],  if  [G]  is  non-singular,  using  classical  routines  (for 
example,  LU  decomposition  and  forward-backward  substitution  are  used  in  the 
present  calculation). 
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Various  choices  of  basis  and  weighting  functions  are  possible.  When  the  two 
functions  are  identical,  wm  =  gm  ,  the  procedure  is  known  as  Galerkin’s  method.  In  the 

present  code,  the  weighting  function  is  different  from  the  basic  function.  The 
weighting  functions  are  a  set  of  delta  functions: 

wm(r)=S(r-rm)  (12) 

where  {rm}  is  a  set  of  points  on  the  conducting  surface.  The  result  is  that  the  integral 
equations  are  projected  onto  Dirac  delta  functions  at  the  points  of  the  discretized 
surface.  The  method  is  referred  as  the  collocation  or  point  matching  method. 

The  basis  functions  chosen  to  expand  the  surface  current  on  each  patch  is  a  set 
of  pulse  functions.  This  selection  simplifies  the  evaluation  of  the  inner  product 
integral  and  ensures  that  the  matrix  [G]  is  well  conditioned.  The  pulse  function 
expansion  for  a  number  of  P  patches  is: 

J5  (r)  =  Z  [A/S  (r/ )+  J 2je2j  (rj  )k  (r)  (13) 

7=1 

where  rj  is  the  point  matching  at  the  centre  of  patch  number  j,  and  Vj(r)  =  1  for  r  on 
patch  j  and  0  otherwise.  The  current  expansion  substituted  into  (7)  leads  to  a  linear  set 
of  equations  that  can  be  written  in  the  matrix  form  (11)  from  where  the  unknown 
expansion  coefficients  Jq  and  hj  representing  the  average  current  density  on  the 
surface  of  the  patch  are  determined. 

The  matrix  element  Gmn  represents  the  tangential  magnetic  field  component  at 
patch  m  due  to  a  surface  current  pulse  on  patch  n.  The  pulse  is  in  the  direction  ei  when 
n  is  odd,  and  in  the  direction  e2  when  n  is  even.  For  m  =  n,  the  contribution  to  the 
vector  product  on  the  flat  surface  is  given  only  by  the  images  of  the  patch  current 
created  on  the  boundaries,  which  is  added  to  the  coefficient  of  Js(r)  in  equations  (7) 
that  is  ±1/2.  In  a  homogeneous  space  the  surface  integral  is  zero  since  the  vector 
product  is  zero.  However,  Gmm  =  ±1/2  in  this  case. 
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4.  Surface  modeling 


The  surface  of  a  perfectly  conducting  (PEC)  object  is  modelled  by  means  of 
multiple  flat  surface  patches  chosen  in  a  manner  to  cover  completely  the  surface  to  be 
modelled.  Thus  the  surface  is  subdivided  into  a  reasonably  number  of  patches 
conforming  as  closely  as  possible  to  both  the  flat  and  curved  surfaces.  The  modelled 
surface  must  be  closed  since  the  patches  only  model  the  side  of  the  surface  from 
which  the  nonnal  is  directed  outward.  The  shape  of  an  elementary  patch  does  not 
affect  the  solution  since  there  is  no  integration  over  the  patch.  The  program  computes 
the  surface  current  in  the  centre  of  each  patch  along  the  orthogonal  unit  vectors  e^r) 
and  e,  (r) ,  which  are  tangent  to  the  surface.  These  current  components  are  essentially 
electric  dipoles  that  radiate  in  the  presence  of  the  interfaces  of  the  stratified  structure 
(4). 

Although  an  elementary  patch  may  be  defined  using  a  shape  (square, 
rectangular,  triangle)  on  the  input  of  the  program,  the  actual  parameters  defining  a 
surface  patch  are  the  Cartesian  coordinates  of  the  patch  centre,  the  components  of  the 
outward  directed  unit  normal  vector  and  the  patch  area. 

As  mentioned  before,  the  particularities  of  the  integral  equation  method  make 
a  model  useable  within  a  frequency  band.  Patch  size  measured  in  wavelengths  is  very 
important  for  accuracy  of  the  results.  A  minimum  of  about  25  patches  should  be  used 
per  square  wavelength  of  surface  area,  with  the  maximum  size  for  an  individual  patch 
about  0.04  square  wavelengths. 

The  optimum  model  for  a  structure  can  be  estimated  by  varying  the  patch 
density  and  observing  the  effect  on  the  results.  Some  dependence  of  results  on 
segmentation  will  always  be  found.  A  large  dependence,  however,  would  indicate  that 
the  solution  has  not  converged  and  more  patches  should  be  used. 
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5.  Numerical  examples 


To  simulate  the  electromagnetic  response  of  a  3-D  PEC  body  in  a  layered 
media,  the  algorithm  based  on  the  Method  of  Moments  integral  equation  was 
implemented  in  a  numerical  code  written  in  Fortran.  The  code  was  tested  against  the 
NEC  code  for  the  particular  case  of  the  free  space  (the  only  capability  of  the  NEC 
code).  However,  in  the  present  code  the  dyadic  Green’s  functions  were  calculated  for 
a  stratified  media  where  the  boundaries  were  placed  far  away  from  the  source,  sensors 
and  scattering  object  to  minimize  their  influence.  The  results  obtained  from  the 
present  code  and  from  the  NEC  code  were  identical  (errors  in  the  fifth  digit). 

The  validity  of  the  program  was  also  checked  by  comparing  results  obtained 
for  the  case  of  an  infinite  long  buried  horizontal  cylinder  excited  by  a  horizontal  line 
source  of  current  above  the  surface  of  a  uniform  half-space  and  parallel  to  the  axis  of 
the  buried  cylinder.  Ogunade  [13]  analyzed  this  2-D  problem  and  obtained  a 
numerical  solution  from  an  extension  of  Dyakonov  [14]  analytical  solution.  The 
approximate  solution  offered  by  Ogunade  [13]  was  obtained  by  truncating  the 
analytical  solution  to  a  finite  number  of  terms. 


Figure  2.  Geometry  of  the  2-D  buried  cylinder  test  problem  (Ogunade  [13]). 
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We  refer  to  Figure  2  for  the  geometry  of  the  test  example.  An  infinite  long 
circular  cylinder  with  radius  ‘a’  and  conductivity  02  is  embedded  in  a  half-space 
conducting  media  having  the  conductivity  ci.  The  free  space  has  zero  conductivity. 
The  cylinder  is  oriented  with  its  axis  along  the  X-axis.  The  calculation  was  made  for 
the  symmetric  case  where  the  infinite  current  carrying  wire  is  located  at  yL  =  zl  =  -h/2. 
The  measurable  quantities  are  the  amplitude  and  phase  of  the  ratio  of  horizontal  to 
vertical  magnetic  fields,  Hy/Hz,  which  are  presented  in  figures  3a  and  3b.  The  figures 

are  defined  in  terms  of  numerical  distances  H  =  h^jcojul)a]  ,  A  =  a^coju(tal  ,  Y  =  y/h 

and  conductivity  contrast  K  =  G2/G1.  In  figures  3a  and  3b  the  numerical  parameters  are 
H  =  0.1,  A  =  0.01,  and  K  =  107.  It  was  considered  that  the  values  of  d  =  lS/m  and  G2 
=  107  S/m  will  offer  a  good  approximation  for  a  perfectly  conducting  cylinder  placed 
in  a  lossy  media.  The  agreement  between  the  two  solutions  is  good  considering  that 
only  5  terms  were  used  in  the  truncated  analytical  solution  and  that  the  coupling 
between  modes  is  neglected  (zero  non-diagonal  tenns). 

An  attempt  to  reproduce  the  scattering  results  presented  in  the  Swedish  FOA 
report  [15]  for  a  PEC  object  in  a  layered  media  generated  differences  of  about  15%, 
which  could  be  explained  by  the  differences  in  solving  for  the  surface  currents  in  the 
two  codes. 
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Figure  3.  Numerical  and  analytical  solutions  of  the  2-D  buried  cylinder. 


The  scattering  calculation  was  applied  to  investigate  the  possibility  to  detect  a 
submerged  submarine  illuminated  by  an  electric  source  placed  on  the  sea  bottom.  The 
source,  which  is  a  100m  long  insulated  wire  with  bare  ends,  carries  10A  sine  wave 
current  with  a  constant  frequency  of  25Hz.  The  receiver  is  a  100m  long  electric  field 
sensor  also  placed  on  the  sea  bottom. 
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The  submarine  is  modeled  with  a  horizontal  cylinder  50m  long,  5m  in 
diameter,  and  has  a  tower  attached.  It  is  submerged  at  10m  depth,  which  represent  the 
distance  between  the  central  axis  of  the  (cylinder)  submarine  and  the  water  level.  The 
surface  of  the  submarine  was  modeled  with  678  patches,  as  illustrated  in  Figure  4 
where  half  of  the  surface  is  presented.  By  increasing  the  number  of  patches  to  1278, 
the  solution  changes  in  its  fourth  digit,  but  the  computation  time  doubled  from  1  to  2 
hours  on  a  2.2  GHz  Pentium  4  computer.  Integration  of  the  tensor  Green’s  functions 
for  the  matrix  elements  and  the  scattering  fields  is  the  most  time-consuming  part  of 
our  routine. 


Figure  4.  Half  of  the  submarine  surface  modeled  with  678  PEC  patches. 

Two  geometries  were  considered,  as  illustrated  in  figures  5a  and  5b:  the 
submarine  is  placed  parallel  with  both  the  source  and  sensor  lines,  but  the  position  of 
the  sensor  relative  to  the  source  is  different.  In  the  first  geometry  (5a),  the  source  and 
the  sensor  lines  are  parallel  at  100m  separation  distance.  In  the  second  geometry,  the 
source  and  the  sensor  are  placed  on  the  same  line  and  they  are  also  100m  separated. 
For  both  geometries  the  water  depth  was  50m  with  a  conductivity  of  3S/m.  The 
sediment  layer  has  the  height  and  conductivity  of  2m  and  0.03S/m  respectively,  and 
the  infinite  rock  bottom  has  the  conductivity  of  0.001  S/m. 
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Electric  field  (uV/m)  Electric  field  (uV/m) 


5a. 


5b. 


Sensor  line  (m)  Sensor  line  (m) 


Figure  5.  The  position  of  the  submarine  in  two  geometries  and  the  scattering  electric  field. 
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Figures  5c  to  5f  presents  the  scattered  electric  field  from  the  submarine  along 
the  sensor  line  for  the  two  geometries.  The  measured  voltage  is  calculated  by 
integration  along  the  sensor  line  and  the  results  are  shown  in  Table  1.  In  this  table  are 
also  presented  the  amplitude  and  phase  of  the  voltage  produced  by  the  direct  field  (in 
the  absence  of  the  submarine). 

The  idea  of  the  active  detection  experiment  is  to  differentiate  between  the 
signal  measured  in  the  absence  of  the  target  (incident  or  primary  fields)  and  the  signal 
with  the  target  present  (primary  plus  secondary  fields).  The  magnitude  of  the 
secondary  field  is  in  general  very  small  in  comparison  to  the  magnitude  of  the  primary 
field  as  it  is  shown  in  Table  1.  For  this  reason,  one  tries  to  maximize  the  ratio  of 
secondary  to  primary  fields  by  choosing  an  appropriate  geometry  for  experiment. 


Table  1.  Direct  and  scattered  voltages  obtained  on  1 00m  sensor  line  at  25Hz 


GEOMETRY 

Vdirect  (fiV) 

^scattered  (|^V) 

RATIO  (%) 

Geo  #1,  D  =  10m 

8.878e+3  Z-176.30 

219.15  Z17.5° 

2.47 

Geo  #1,  D  =  50m 

8.878e+3  Z-176.30 

253.09  Z14.4° 

2.85 

Geo  #1,  D  =  90m 

8.878e+3  Z-176.30 

217.35  Z17.3° 

2.44 

Geo  #2,  D  =  0m 

1.473e+3  Z47.20 

27.27  Z-116.10 

1.85 

Geo  #2,  D  =  100m 

1.473e+3  Z47.20 

23.83  Z  -90.6° 

1.61 

Geo  #2,  D  =  200m 

1.473e+3  Z47.20 

26.74  Z-116.40 

1.81 
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6.  Conclusions 


A  numerical  method  was  developed  for  determining  the  low  frequency 
response  of  perfectly  conducting  (PEC)  submerged  or  buried  3-D  objects  due  to  an 
electric  or  magnetic  source.  Both  the  scattering  object  and  the  source  are  located  in  a 
layered  media.  The  surface  of  the  object  is  partitioned  into  small  flat  patches.  Each 
patch  is  defined  by  its  centre  and  the  surface  normal  with  respect  to  a  known 
coordinate  system.  The  scattering  problem  is  formulated  in  terms  of  the  magnetic  field 
integral  equation  (MFIE).  The  magnetic  Green’s  function  appearing  in  this  equation 
was  calculated  for  a  layered  medium.  Then  the  Method  of  Moments  is  used  to 
transfonn  the  integral  equation  into  a  matrix  equation.  We  wish  to  point  out  that  the 
matrix  fill  time,  not  the  available  computer  memory,  is  the  overriding  factor  that  puts 
practical  limits  on  the  size  of  objects  that  can  be  analysed  by  the  technique  presented 
here. 

The  total  measurable  electric  field  was  separated  into  two  terms:  (a)  the  direct 
(primary)  field,  which  is  the  field  as  if  the  scattering  object  were  absent,  and  (b)  the 
anomalous  (secondary)  field,  which  represented  the  presence  of  the  object.  The 
calculations  we  have  shown,  in  addition  to  verifying  our  numerical  solution,  indicate 
that  object  response  is  strongly  influenced  by  the  presence  of  layering.  Results  were 
presented  for  a  target  located  in  a  conducting  layer  of  seawater.  The  scattered  EM 
field  is  strongly  dependent  on  the  orientation  of  the  vessel  with  respect  to  the  source 
and  receiver.  However,  the  scattered  intensity  is  typically  a  few  percent  of  the  direct 
field.  For  a  good  detection  one  needs  to  maximize  the  ratio  of  secondary  field  to 
secondary  plus  primary  fields,  which  could  be  achieved  by  an  appropriate  setting  of 
experiment. 


16 


DRDC  Atlantic  TM  2004-246 


7.  References 


[1]  W.  L.  Ko  and  R.  Mittra,  Extremely  low  frequency  modelling  in  lossy  media  using 
FDTD  with  application  in  seafloor  characterisation,  Electromagnetics,  vol.  15,  p.  587- 
602, 1995. 

[2]  P.  E.  Wannamaker,  G.  W.  Hohmann  and  W.  A.  SanFilipo,  Electromagnetic 
modeling  of  three  dimensional  bodies  in  layered  earths  using  integral  equations, 
Geophysics,  vol.  49  (1),  p.  60-74,  1984. 

[3]  K.  A.  Michalski  and  J.  R.  Mosig,  Multi-layered  media  Green's  functions  in  integral 
equation  fonnulations,  IEEE  Trans,  on  antennas  and  propagation,  vol.  45  (3),  p.  508- 
519,  1997. 

[4]  T.  Jun  Cui,  W.  Wiesbeck  and  A.  Herschlein,  Electromagnetic  scattering  by 
multiple  three-dimensional  scatterers  buried  under  multilayered  media  -  Part  I: 
Theory,  IEEE  Trans,  on  Geoscience  and  Remote  Sensing,  vol.  36  (2),  p.  526-534, 
1998. 

[5]  T.  Jun  Cui,  W.  Wiesbeck  and  A.  Herschlein,  Electromagnetic  scattering  by 
multiple  three-dimensional  scatterers  buried  under  multilayered  media  -  Part  II: 
Numerical  implementations  and  results,  IEEE  Trans,  on  Geoscience  and  Remote 
Sensing,  vol.  36  (2),  p.  535-546,  1998. 

[6]  N.  Geng,  C.  E.  Baum  and  L.  Carin,  On  the  low-frequency  natural  response  of 
conducting  and  penneable  targets,  IEEE  Trans,  on  Geoscience  and  Remote  Sensing, 
vol.  37(1),  p.347-359,  1999. 

[7]  S.  S.  Bindiganavale  and  J.  L.  Volakis,  Comparison  of  three  FMM  techniques  for 
solving  hybrid  FE-BI  systems,  IEEE  antennas  and  propagation  magazine,  vol.  39  (4), 
p.  47-50,  1997. 

[8]  X.  Q.  Sheng,  J.  M.  Jin,  J.  Song,  C.  C.  Cheng  and  W.  C.  Chew,  On  the  formulation 
of  hybrid  finite-element  and  boundary-integral  methods  for  3-D  scattering,  IEEE 
Trans,  on  antennas  and  propagation,  vol.  46  (3),  p.  303-3 11,  1998. 

[9]  R.  C.  Ireland,  P.  C.  Coles,  D.  Rodger  and  P.  J.  Leonard,  Low  frequency 
electromagnetic  scattering  by  a  submerged  finite  length  cylinder,  Proceedings  of 
International  Conference  on  Marine  Electromagnetics,  Marelec  97,  Brest-France,  June 
1997. 


DRDC  Atlantic  TM  2004-246 


17 


[10]  A.  Sahin  and  E.  L.  Miller,  Recursive  T-matrix  methods  for  scattering  from 
multiple  dielectric  and  metallic  objects,  IEEE  Transactions  on  Antennas  and 
Propagation,  vol.  46  (5),  p.  672-678,  1998. 

[11]  M.  Birsan,  "VLF  code  for  horizontally  stratified  media,  Defence  Research 
Establishment  Atlantic  Technical  Memorandum,  DREA  TM  2000-151,  December 
2000. 

[12]  J.  T.  Weaver,  “The  quasi-static  field  of  an  electric  dipole  embedded  in  a  two- 
layer  conducting  half-space”,  Canadian  Journal  of  Physics,  vol.  45,  pp.  1981-2002, 
1967. 

[13]  S.O.  Ogunade,  Electromagnetic  response  of  an  embedded  cylinder  for  line 
current  excitation,  Geophysics,  vol.  46  (1),  p.  45-52,  1981. 

[14]  B.  P.  Dyakonov,  The  diffraction  of  electromagnetic  waves  by  a  circular  cylinder 
in  a  homogeneous  half-space,  Izvestiya,  Acad.  Sci.  USSR,  Geophysics  ser.  No.  9,  p. 
950-955,  1959. 

[15]  J.  Mattsson,  A  computer  program  for  EM-scattering  by  smooth  3D  objects  in 
horizontally  stratified  media,  Swedish  Defence  Research  Establishment  technical 
report  FOA-99-409,  Stockholm,  Sweden,  1999. 


18 


DRDC  Atlantic  TM  2004-246 


List  of 

symbols/abbreviations/acronyms/initialisms 


DND 

EM 

PEC 


Department  of  National  Defence 

Electromagnetic 

Perfectly  conducting 


DRDC  Atlantic  TM  2004-246 


19 


Distribution  list 


Library  DRDC  Atlantic  -  6  copies 
Marius  Birsan  -  1  copy 
DRDKIM  -  1  copy 


20 


DRDC  Atlantic  TM  2004-246 


DOCUMENT  CONTROL  DATA 

(Security  classification  of  title,  body  of  abstract  and  indexing  annotation  must  be  entered  when  the  overall  document  is  classified) 

1.  ORIGINATOR  (the  name  and  address  of  the  organization  preparing  the  document.  2.  SECURITY  CLASSIFICATION 

Organizations  for  whom  the  document  was  prepared,  e.g.  Centre  sponsoring  a  (overall  security  classification  of  the  document 

contractor's  report,  or  tasking  agency,  are  entered  in  section  8.)  including  special  warning  terms  if  applicable). 

DRDC  Atlantic  UNCLASSIFIED 

9  Grove  Street 

Dartmouth,  NS,  Canada  B2Y  3Z7 

3.  TITLE  (the  complete  document  title  as  indicated  on  the  title  page.  Its  classification  should  be  indicated  by  the  appropriate 
abbreviation  (S,C,R  or  U)  in  parentheses  after  the  title). 

Modeling  the  electromagnetic  scattering  by  3-D  objects  in  a  layered  media  (U) 


4.  AUTHORS  (Last  name,  first  name,  middle  initial.  If  military,  show  rank,  e.g.  Doe,  Maj.  John  E.) 

Marius  Birsan 

5.  DATE  OF  PUBLICATION  (month  and  year  of  publication  of 
document) 

December  2004 

6a.  NO.  OF  PAGES  (total 

containing  information  Include 
Annexes,  Appendices,  etc). 

17 

6b.  NO.  OF  REFS  (total  cited 
in  document) 

15 

7.  DESCRIPTIVE  NOTES  (the  category  of  the  document,  e.g.  technical  report,  technical  note  or  memorandum.  If  appropriate,  enter  the 
type  of  report,  e.g.  interim,  progress,  summary,  annual  or  final.  Give  the  inclusive  dates  when  a  specific  reporting  period  is  covered). 

TECHNICAL  MEMORANDUM 

8.  SPONSORING  ACTIVITY  (the  name  of  the  department  project  office  i 

Defence  R&D  Canada  -  Atlantic 

PO  Box  1012 

Dartmouth,  NS,  Canada  B2Y  3Z7 

or  laboratory  sponsoring  the  research  and  development.  Include  address). 

9a.  PROJECT  OR  GRANT  NO.  (if  appropriate,  the  applicable  research  9b.  CONTRACT  NO.  (if  appropriate,  the  applicable  number  under 
and  development  project  or  grant  number  under  which  the  document  which  the  document  was  written), 

was  written.  Please  specify  whether  project  or  grant). 


Project  1 1  cj 

10a  ORIGINATOR'S  DOCUMENT  NUMBER  (the  official  document  10b  OTHER  DOCUMENT  NOs.  (Any  other  numbers  which  may  be 
number  by  which  the  document  is  identified  by  the  originating  assigned  this  document  either  by  the  originator  or  by  the 

activity.  This  number  must  be  unique  to  this  document.)  sponsor.) 

DRDC  Atlantic  TM  2004-246 

11.  document  availability  (any  limitations  on  further  dissemination  of  the  document,  other  than  those  imposed 
by  security  classification) 

(  X  )  Unlimited  distribution 

(  )  Defence  departments  and  defence  contractors;  further  distribution  only  as  approved 

(  )  Defence  departments  and  Canadian  defence  contractors;  further  distribution  only  as  approved 

(  )  Government  departments  and  agencies;  further  distribution  only  as  approved 

(  )  Defence  departments;  further  distribution  only  as  approved 

(  )  Other  (please  specify):  “Controlled  by  Source” 

12.  DOCUMENT  ANNOUNCEMENT  (any  limitation  to  the  bibliographic  announcement  of  this  document.  This  will  normally  correspond  to  the 
Document  Availability  (11).  However,  where  further  distribution  (beyond  the  audience  specified  in  (11)  is  possible,  a  wider  announcement 
audience  may  be  selected). 


DRDC  Atlantic  mod.  May  02 


13.  ABSTRACT  (a  brief  and  factual  summary  of  the  document.  It  may  also  appear  elsewhere  in  the  body  of  the  document  itself.  It 
is  highly  desirable  that  the  abstract  of  classified  documents  be  unclassified.  Each  paragraph  of  the  abstract  shall  begin  with  an 
indication  of  the  security  classification  of  the  information  in  the  paragraph  (unless  the  document  itself  is  unclassified)  represented 
as  (S),  (C),  (R),  or  (U).  It  is  not  necessary  to  include  here  abstracts  in  both  official  languages  unless  the  text  is  bilingual). 

The  underwater  electromagnetic  (EM)  active  detection  technique  attempts  to  reveal  a 
stealth  iron-hull  vessel  (submarine  or  surface  warship)  located  in  a  horizontal  layer  of 
seawater.  Our  goal  is  to  numerically  investigate  the  EM  scattering  process  in  the  presence 
of  a  metallic  vessel  of  arbitrary  shape  placed  inside  a  lossy  media.  To  solve  the  EM 
scattering  problem  involving  perfectly  conducting  (PEC)  objects  we  make  use  of  the 
Method  of  Moments  (MoM)  integral  equation  technique.  The  EM  scattering  numerical 
code  was  written  in  Fortran.  The  solution  is  applied  to  the  investigation  of  EM  scattering 
by  a  submarine  modelled  as  a  tapered  long  cylinder  plus  the  tower  when  is  placed  in 
different  positions  relative  to  the  source  and  the  sensors.  It  is  shown  that,  by  an 
appropriate  design  of  the  experiment,  the  signal  indicating  the  presence  of  the  target  can 
be  increased  relative  to  the  signal  obtained  in  its  absence. 


14.  KEYWORDS,  DESCRIPTORS  or  IDENTIFIERS  (technically  meaningful  terms  or  short  phrases  that  characterize  a 
document  and  could  be  helpful  in  cataloguing  the  document.  They  should  be  selected  so  that  no  security  classification  is 
required.  Identifiers,  such  as  equipment  model  designation,  trade  name,  military  project  code  name,  geographic  location  may 
also  be  included.  If  possible  keywords  should  be  selected  from  a  published  thesaurus,  e.g.  Thesaurus  of  Engineering  and 
Scientific  Terms  (TEST)  and  that  thesaurus-identified.  If  it  not  possible  to  select  indexing  terms  which  are  Unclassified,  the 
classification  of  each  should  be  indicated  as  with  the  title). 


electromagnetic  scattering,  layered  lossy  media,  active  vessel  detection,  method  of  moments 


DRDC  Atlantic  mod.  May  02 


This  page  intentionally  left  blank. 


Defence  R&D  Canada  R  &  D  pour  la  defense  Canada 


Canada's  leader  in  defence 
and  National  Security 
Science  and  Technology 


Chef  de  file  au  Canada  en  matiere 
de  science  et  de  technologie  pour 
la  defense  et  la  securite  nationale 


DEFENCE 


www.drdc-rddc.gc.ca 


